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Baryogenesis with Higher Dimension Operators 
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We propose a simple model of baryogenesis comprised of the standard model coupled to a singlet X 
via higher dimension operators O. In the early universe, X is thermalized by O mediated scattering 
processes before it decouples relativistically and evolves into a sizable fraction of the total energy 
density. Eventually, X decays via O in an out of equilibrium, baryon number and CP violating 
process that releases entropy and achieves baryogenesis for a broad range of parameters. The decay 
can also produce a primordial abundance of dark matter. Because X may be as light as a TeV, 
viable regions of parameter space lie within reach of experimental probes of n-n oscillation, flavor 
physics, and proton decay. 



I. Introduction 

The standard model (SM) cannot explain the observed 
matter-antimatter asymmetry of the universe, and so new 
physics is required. In this letter we propose a simple 
scenario for baryogenesis consisting of the SM plus an 
inert multiplet of states X . These states interact weakly 
with the SM through baryon number and CP violating 
higher dimension operators O set by the scale A. 

The process of baryogenesis occurs in the four stages 
depicted in Fig. [TJ In the beginning, 

i) X is thermalized with the SM plasma. 

This condition is possible provided Tr, the reheating 
temperature, is greater than mx, the mass of X. Hence, 
thermalization occurs automatically via scattering in the 
SM plasma mediated by O or the ultraviolet dynamics 
which generates O. Once the universe cools sufficiently, 
O mediated scattering goes out of equilibrium and 

ii) X decouples relativistically from the SM plasma. 

Once X leaves equilibrium, it redshifts like radiation until 
temperatures drop below mx , at which point X becomes 
non-relativistic. Once X begins to redshift like matter, 

Hi) X evolves into a large fraction of the total energy. 

During this period the energy density in X is greater 
than that of any given relativistic species, and may even 
come to dominate the total energy density, sending the 
universe into a matter dominated phase. The epoch of 
X domination terminates when 

iv) X decays, yielding a primordial baryon asymmetry. 

Crucially, these out of equilibrium decays of X occur via 
the very same baryon number and CP violating higher 
dimension operators O that initially thermalize X in the 
early universe. Interference between tree and one-loop 
decay amplitudes generate a baryon asymmetry in the 
final state, as depicted in Fig. [5] for an explicit model. In 
certain models, X decays can also generate a primordial 
abundance of dark matter (DM). 

Let us highlight the key features of this baryogenesis 
scenario. First, since this mechanism allows for low scale 



baryogenesis, the operators O may be indirectly probed 
through n-n oscillations, flavor violation, and proton de- 
cay. Second, this setup is quite minimal, in that only 
a handful of new particles X are required and the very 
same operators O that produce X initially also mediate 
its decay. As we will see later, a subset of the X particles 
can even be DM. Third, this setup exploits a cosmological 
"fixed point" arising because X is typically thermalized 
for a very broad range of reheating temperatures. 

As is well-known, the production and thermalization 
of inert particles at reheating is a ubiquitous difficulty in 
theories beyond the SM. This issue arises in the cosmol- 
ogy of gravitinos QHp , axinos [3-H^. photini [fj, [l4| . 
and goldstini [la . Il6j . Transforming this peril into a 
blessing is an old idea, e.g. in models linking gravitino 
or axino domination to baryogenesis in R-parity violat- 
ing supersymmetry |17l.ll8|. However, we argue that this 
mechanism applies much more broadly and is a natu- 
ral byproduct of additional singlet states coupled to the 
SM via baryon number and CP violating higher dimen- 
sion operators — the out of equilibrium condition arises 
from relativistic decoupling and decays of X. Alterna- 
tively, the out of equilibrium condition for X can be 
achieved through heavy particle decays |19M23| . first or- 
der phase transitions |24l427| . or rolling scalars [2a-l3fj||. 
Mechanisms involving higher dimension operators have 
also been discussed in more specific contexts [31J, [32j ■ 

In Sec. II we present a simple example model that illus- 
trates the salient features of our setup. We then discuss 
constraints from experimental limits in Sec. Ill and con- 
clude in Sec. IV. 

II. The Model 

In this section we present a simple theory illustrating 
our mechanism for baryogenesis. Consider the SM aug- 
mented by a multiplet of gauge singlet Majorana fermions 
Xj with mass m/. The interaction Lagrangian for Xj is 



~X r 



{X I u i )(Xju j ) 



A 



li'jk 

A 2 



(X lUi )(djd k ) + c.c, (I) 



where i,j, k label right-handed quark flavors. Lorentz 
indices are contracted implicitly among terms in paren- 
theses, while color indices are contracted implicitly in the 
unique way. It is straightforward to include other Lorentz 
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FIG. 1: The four stages of baryogenesis, shown in terms of 
the evolution of the energy density in SM radiation and X as 
a function of scale factor. The decay of X may occur before 
or after X grows to dominate the total energy density. 



and flavor structures into the Lagrangian, but such terms 
will not qualitatively alter the mechanics of the model. 

In terms of symmetries, baryon number is violated be- 
cause Xi are Majorana, while CP is violated because the 
couplings Kjjij and Xnjk are complex. Note that there is 
an exact, unbroken Z2 subgroup of baryon number under 
which Xj and the quarks are all odd. 

Let us now describe the cosmological history of this 
model. To begin, we assume that the SM is reheated to a 
temperature Tr > mj shortly after inflation. If Tr > A, 
then the effective theory described in Eq. ((1) does not 
apply, but any renormalizable ultraviolet completion of 
these higher dimension operators will generically induce 
tree level scattering processes that thermalize X. On 
the other hand, if Tr < A, then the higher dimen- 
sion operator description is valid, and the interactions 
in Eq. (Q]) will mediate high energy scattering processes 
such as Uiflj —} XiXj, Uidj —} Xidk, djdk — > Xiu% which 
also tend to thermalize Xj . The thermally averaged pro- 
duction cross-section for Xj scales as ((tv)j ~ cjT 2 /A 4 , 
where the proportionality factor cj depends on Xnjk and 
kjj^. Thus, X scattering is dominated by ultraviolet 
processes, and is most important at Tr. This effect is 
familiar from supersymmetric cosmology, where overpro- 
duction of gravitinos during reheating places a stringent 
limit on Tr [l|. Similar limits have been computed for a 
general hidden sector cosmology [33j . 

The critical decoupling temperature T/j 7 defines the 
temperature at which these scattering processes go out 
of equilibrium, i.e. when n cq (av)i ~ H where n Gq is the 
equilibrium number density of Xj and H is the Hub- 
ble parameter. Together with the scaling of {crv)i, this 
implies that Td, ~ (A 4 /mpi) 1 / 3 where mpi ~ 2.4 x 
10 18 GeV is the reduced Planck mass. In summary, if 
T R > A or A > Tr > T Dl , then Xi will be thermalized 
during reheating. This thermalization condition is eas- 
ily satisfied for sufficiently high values of Tr, which we 
assume for the remainder of our discussion. 



While Xj is thermalized initially, it leaves equilibrium 
once temperatures drop below Tjj 1 . After decoupling, 



the yield of Xi is given is given by 

Yi - =#^, (2) 

s{T Dl ) 

where s is the entropy density. The yield is constant 
in the absence of entropy production. Because we are 
interested in a case in which Xj decouples while it is 
relativistic, we assume throughout that T]j I > mj. 

Once temperatures drop below mi, Xi becomes non- 
relativistic and its energy density begins to redshift like 
matter, since pi(T)/s(T) — mjYj is a constant. From 
this point forward, the energy density of X\ will evolve 
to dominate that of any relativistic species. During this 
era, Xi may even come to dominate the total energy 
density, at which point the universe will enter a matter 
dominated phase. 

This period of Xi domination ends when Xj decays 
via processes of the form Xj — >• Ujdjdk, Ujdjdk, XjUjUj. 
This final state of baryogenesis is similar to that of [23| . 
The associated partial decay widths are 



r(Xj 



Uidjdk) 



r(X 7 -> Xju iUj ) 



\Xnjk\ 2 m 5 j 
512tt 3 A 4 ' 
\Kijij\ 2 m 5 j 



(3) 

(4) 



1024tt 3 A 4 ' 

ignoring kinematic factors arising from masses of final 
state particles. The lifetime of Xj is constrained by 
number of cosmological constraints. First, if A is too 
high, the model is constrained by stringent limits from 
big bang nucleosynthesis (BBN) |34l . l35| on late time in- 
jection of electromagnetic energy. The lifetime of Xi is 
thus bounded by 77 < 1 s, where 
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(5) 



and we have defined the effective couplings Xj = 
J2ijk l^Ujk \ 2 +J2jij \ K Uij I 2 / 4 where the sums range over 
kinematically allowed final states. Eq. (JSJ) demonstrates 
that BBN bounds are satisfied for a broad range of pa- 
rameter space. Second, if A is too low, then (av)i will be 
large and scattering may keep Xj in equilibrium down 
to temperatures of order to/. In this case, Tjj [ < m.j 
and Xj decouples non-relativistically. While a residual 
baryon asymmetry maybe still persist, a correct evalua- 
tion would require a full analysis of Boltzmann equations 
which goes beyond the scope of this work, so we restrict 
to the case where Xj decouples relativistically. 

All but the lightest of the Xj will have CP violating 
decay modes of different baryon number, so their decays 
produce a final state baryon asymmetry through one-loop 
interference. The asymmetric width of Xi — > uidjd^ 
is given by interference between tree and loop diagrams 
depicted in Fig. [5] Ignoring kinematic factors in the final 
and intermediate states, 



r(A/ ->• Uidjdk) - r(Xi ->■ u t djd k ) 
Im(AJy fe KjjKAjy fc ) m\ 



E 



51207T 4 



A6 



(6) 



where here the sums range over kinematically accessible 
final and intermediate states. We define an asymmetry 
parameter for each Xj decay by 



a = J2 B(f)\BR(Xi "> /) - BRPG ->- /")] 



/ 

1 Si m 2 
20^ Af A 1 "' 



(7) 



where / sums over final states, BR denotes the branching 
ratio of a given process, and B is the baryon number 
of each final state. Here we have defined the quantity 
5i = J2jijkiI m (^*iijk K uii^Jijk) to characterize the net 
CP violation associated with A/. 

In order to compute the net baryon asymmetry, let us 
first consider the decay of a single component, Xj. The 
cosmology depends sensitively on the relative size of pj, 
the energy density in Xj, and pr, the total energy den- 
sity in radiation, evaluated just prior to decay. If Xi 
decays very soon after it becomes non-relativistic, then 
its energy density is of order that of a single relativis- 
tic species, which we dub the "weak domination regime", 
Pi *C pr. Little entropy is produced and the temperature 
of the radiation remains more or less constant. However, 
if Xj decays quite late, then it dominates the total energy 
density of the universe, which we dub the "strong domi- 
nation regime", pi 3> pr. Thus, A/ decays will boost the 
temperature of the radiation bath to an effective temper- 
ature Ti determined by the total energy density injected, 
pi = Tr 2 g*Tf/30 = 3ff 2 m P1 when ti ~ 1/H, so 



T, 



90 



ir 2 g*(Ti) 



1/4 



a/topi/t/, 



(8) 



where g* counts relativistic degrees of freedom. 

The asymmetric baryon number generated by Xj de- 
cays is given by 



?]i = eiYidi, 



(9) 



where Yj is defined in Eq. (J5|) and the dilution factor di 
is the ratio of the entropy density before and after Xi 
decays, so 



1 , Pi < PR 

3& . Pi* PR 



(10) 



so the dilution factor is much smaller in the strong dom- 
ination regime. Applying Eqs. |(5J), ^ and (fTU|) we find 
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(11) 



in the weak domination regime, pi <C pr- Meanwhile, 
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FIG. 2: Tree and one-loop diagrams for Xi — > Uidjdk which 
interfere to produce a primordial baryon asymmetry. 



in the strong domination regime, pi 3> pr. Note that 
sphaleron processes will partially wash out the baryon 
asymmetry if Tj > mw/ctw, where mw is the W boson 
mass, a>w = g 2 /4ir, and g is the SU(2)i gauge coupling. 
In this case, the net baryon asymmetry is processed ac- 
cording to m -4 (28/79)77/. 

With the expressions in Eqs. (|TTj) and p^|) . it is 
straightforward to compute the net baryon asymmetry 
generated from all Xi decays. If the masses and cou- 
plings are not hierarchical, then each Xj should decay 
around a similar time. In this case, the Xi should ei- 
ther all be in the weak domination regime or all be in 
the strong domination regime. For the former, little en- 
tropy is produced by each Xi decay, and the net baryon 
asymmetry is simply given by the sum of all r\i . For the 
latter, entropy is substantially produced in each decay, 
thus diluting the asymmetry generated in earlier epochs. 
In this case the net baryon asymmetry is dominated by 
77/ from latest of the Xj decays. 

Finally, let consider the issue of DM. Let us denote a 
stable component of the Xi multiplet by Adm- If Adm 
is lighter than the proton, then Adm is exactly stable 
because it is the lightest odd particle under the unbro- 
ken Z2 subgroup of baryon number. The primordial relic 
abundance of Adm is ^dm = "Idm^dm(s/Vc)oi where 
(p c /s) Q ~ 3.6 h 2 x 10~ 9 GcV for h ~ 0.67 [3|. The DM 
abundance arises from two sources, Ydm = ^dm + -^dm 
arising from thermal scattering during initial reheat- 
ing and decays of heavier A/, respectively. Assuming 
that baryogenesis is dominated by the decays of a single 
species Xi, these contributions are 






m y"pM 

ei Yj ' 



Y^ = ^BR(A 7 



Adm), 



(13) 
(14) 



(12) 



where Ydm and Yi are as defined in Eq. ([2]). Typically, 
Ydm and Yi will be comparable, and BR(A/ -4 Adm) ~ 
O(l), so the contributions to the DM abundance from 
thermal scattering and decays will be of similar order, 
but both ~ 1/e/ larger than the asymmetric yield. 

Finally, we summarize the allowed parameter space for 
a single species of A/ in Fig. [3] The grey region indi- 
cates where the higher dimension operator description is 
invalid because mi > A. For A/ = 1, the blue region 



10' 



10' 



10 



> 10 

CO 

< 10 8 



10 B 



10" 



. BBN limits „ 


1 


dF s * 




> . 


// 


''a 1 ( 






'' 






" j7'\ 
















'/V \^ / 


nvalid 


- 


/ «.& / 








: & / A 


/ = 


Adm 


= 1 






5/ 


= 5 


/ BR(X 7 -> 


X, 


dm) = 


0.5 



10 4 10 6 10 8 10 10 10 12 10 14 



m 7 (GeV) 

FIG. 3: Cosmologically allowed regions of parameter space. 
The grey region lies outside the regime of the effective theory. 
The blue region is disfavored by BBN, while the purple region 
is excluded by the requirements that Xi decouple relativisti- 
cally and washout be evaded. The yellow region accommo- 
dates the observed primordial baryon asymmetry, while the 
green dashed lines denote the required DM mass in keV. 



depicts the parameter space excluded by BBN limits on 
late decays of A/, while the purple region depicts the pa- 
rameter space excluded by requiring that Xj decouples 
relativistically — i.e. it is not thermalized by scattering 
processes at temperatures of order to,/. The purple re- 
gion is very similar to the region excluded by washout 
from scattering processes, assuming c/ = A//47T for all 
Xj. Furthermore, taking that Si = 5 and that Xj is 
the primary origin of baryogenesis, then the yellow band 
indicates where 10~ u < m < 10~ 10 . Note that this 
choice for Si is not in a strong coupling regime because 
all couplings are really normalized to a higher dimension 
operator scale A. As noted earlier, the model also has the 
option of including primordial relic DM. Requiring that 
Udm'i 2 — 0.11 [36| fixes the DM mass, which is denote by 
green dashed lines for todm =0.1 and 1 keV. Thus, for 
0(1) couplings, the observed baryon asymmetry is gener- 
ated in the regime in which our effective theory analysis 
is valid, and DM can also be accommodated. 

III. Experimental Signatures 

Our proposal offers experimentally observable conse- 
quences connected with the operators directly involved 
in asymmetry generation. Baryon number violating op- 
erators will typically induce highly constrained n-n os- 
cillations via the effective operator (ud)(dd)(ud) / M^_ n , 
where u and d are the right-handed up and down quarks, 
while Lorentz indices are contracted within the parenthe- 
ses and color indices are contracted in the unique way. 
This operator is not induced at tree level in the model de- 
fined in Eq. (JT|) , due to an accidental antisymmetry in the 
flavor indices of the coupling constant Xnjk- However, 
this operator will be induced at loop order, and more 
generally will be present at tree level if there are higher 
dimension operators in addition to those in Eq. ([T]). For 
example n-n oscillation will be induced if there are oper- 
ators of the form \' Ii - k (Xidj){dkUi) / A 2 . Integrating out 



X\ will produce the n-n operator with an effective cut- 
off M^_ n ~ A 4 TO//Aj 2 m . The characteristic time scale 
of n-n oscillation goes as r n _ n ~ M^_^/(3 x 10~ 4 GeV 6 ) 
[371 , l38fl . which together with the experimental bound, 
T n -n > 2.4 x 10 s s (3^ implies 



A > 3.2 x 10 6 GeV \X' nil 



,1/2 



!TeV\ 1/4 



mi 



(15) 



so n-n oscillations could offer a sensitive probe of the low 
scale variants of this baryogenesis mechanism. 

Flavor violation offers another possible probe of this 
model. In particular, K°-K° mixing is mediated by 
the operator (dd)(ss)/M?, _^- , which is induced at loop- 



16 7 t 2 A 4 /A7 111 A^ 122 to 2 . Corn- 



level, where Af 2 ^ 

paring the estimated mixing rate with the experimental 

bound, ImA/12 < 3.3 x lfr 18 GeV J4(Jfi3, gives 



A > 4.4 x 10 4 GeV Im(Aj in A^ 



a/4 



(t 



mi 



TcV 



1/2 



,(16) 



which can be competitive with n-n limits. 

In theories where there exists a cosmologically sta- 
ble dark matter candidate Adm, there are stringent 
limits on proton decay via the process p — > 7t + Adm, 
whose decay rate is estimated as T(p —> 7t + Adm) ~ 
•^dm^t-^qcdA^ttA 4 , where m p is the proton mass and 
Aqcd ~ 250 McV is the QCD scale 0|. Then the ex- 
perimental bound, T p ^^+y > 2.5 x 10 31 yr [4fJ|, gives a 
very stringent limit on the cutoff 



A > 5.5 x 10 14 GeV A 



1/2 

DM 



( A 



QCD 



V 250 McV 



(17) 



In order to evade the proton decay bound for the DM 
model, we must assume a hierarchical flavor structure 
in the coupling Aqm to the light quarks. This can 
be accommodated in models of minimal flavor violation 
(MFV) , e.g . in R-parity violating supersymmetric theo- 
ries |38l.[45|. 

IV. Conclusions 

In this letter we propose a simple baryogenesis model 
comprised of the standard model coupled via higher di- 
mension operators to a neutral multiplet X. In the early 
universe, X will typically be thermalized after reheating, 
but it eventually decouples and becomes non-relativistic. 
Subsequently, X evolves into a sizable fraction of the to- 
tal energy density of the universe until it decays out of 
equilibrium, yielding a baryon asymmetry and possibly 
dark matter. This scenario can be experimentally probed 
via n-n oscillations, flavor violation, and proton decay. 
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